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I. INTRODUCTION 


The coupling of a mixed layer model with an ocean acoustic 
model is an opportunity to make acoustic predictions from earlier BT 
information and meteorological analyses and forecasts. The thermo- 
dynamic model is initialized with an observed temperature profile, 
and it makes use of meteorological information that is currently 
available such as wind speed and solar radiation. The thermody- 
namic forecasts are then input into an acoustic model. In the 
past, similar studies have concentrated on one particular acoustic or 
thermal model. R. H. Fisher (1981) investigated the variability and 
sensitivity of a coupled model system. : He found that a one dimen- 
sional thermal model ¡integrated in: time at a single point (Ocean 
Station Papa 50°N 145°W) predicted mixed layer structure better 
than did the Expanded Ocean Thermal Structure (EOTS) system 
which was currently in use at the Fleet Numerical Oceanography 
Center (FNOC). 

The working hypothesis for this research is that it is possible 
to make accurate acoustic forecasts using recent BT information to 
initialize an ocean mixed layer model driven by atmospheric forcing. 
` The atmospheric forcing is readily available from atmospheric anal- 
ysis and predictions by FNOC. The Garwood model is a proven 
mixed layer model which employs a mathematical model for turbulent 
entrainment (Martin, 1985). The results from the Garwood model 
are input into the RAYMODE Passive Propagation Loss Computer 
Program which uses the RAYMODE propagation loss model, producing 
an acoustic forecast. The RAYMODE model uses aspects of both 
ray and normal mode theory to calculate acoustic pressure. Using 
this system, the acoustic forecaster would not be constrained to 
making predictions based on climatological data alone. 

The purpose of this thesis is to test the working hypothesis 


by expanding the originally one dimensional scope of this research 


to include a line of data points, making the acoustic analysis two 
dimensional. Also proposed is a change in the method Fisher used 
to merge upper layer thermodynamic model output with deeper ocean 
historical information. The thermodynamic model being used 
provides an output of temperature every 5 m down to 200 m. The 
original approach (Fisher, 1981) was to connect the final output 
point at 200 m to historical data available at FNOC. In almost 
every case, there was a large and unrealistic discontinuity in the 
resulting thermal profile. This manifested itself in a large sound 
velocity gradient when input into the acoustic model, and this 
anomaly tended to dominate the acoustic model output. 

The data used in this research was acquired in support of the 
joint U.S.-Canadian Storm Transfer and Response Experiment 
(STREX) which was held in the fall of 1980. The general objec- 
tives of STREX were to understand the physical processes of the 
boundary layers of the atmosphere and ocean in mid-latitude storms, 
the interactions of the two boundary layers, and the interactions 
with larger-scale phenomena  (Miyake, 1980). Observations were 
made in the vicinity of Ocean Station PAPA as storms passed 
during the fall of 1980. The acoustic portion of the experiment 
was called ASTREX and was designed to supplement STREX with 
acoustic measurements relating to the study of the effects of winter 
storms (Dunlap and Andersen, 1985). The data used in this 
research consisted of airborne expendable bathythermograph (AXBT) 
information acquired from a series of six P-3 flights between 15 
November and 5 December 1980. These flights were carried out as 
a part of the Naval Postgraduate School's research effort during 
STREK. Each of the nine hour missions flew northwestward from 
Cape Mendocino, California outbound along the center track (Fig. 
E: The spacing between the stations was 55.6 kilometers (30 
nm). 

The complete set of navigation equipment on tle P-3C was used 
in calculating the position of each of the deployed AXBTs and the 


cumulative error of the inertial navigation system was checked and 
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Figure 1.1 ASTREX center track AXBT stations. 


recorded. This error was less than 4 nautical miles for the flights 
selected for this research (Lundell, 1981). This data is very suit- 
able for thermal-acoustic model analysis because it allows for 


comparison between thermal model output and actual “AXBT 
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observations made at the same time. Also, because the stations are 
in a straight line and are equally spaced, the transition from one 


dimension to two dimensions is easily accomplished. 


II. OCEANIC CHARACTERISTICS OF THE PROJECT AREA 





The region chosen for the Naval Postgraduate School's portion 
of the STREX project was that portion of the northeast Pacific 
Ocean between latitudes 40°N and 50°N and 126°W and 139°W. The 
area is located between the Pacific Subarctic Water Mass and the 
Pacific Equatorial Water Mass (Tabata, 1965). The North Pacific 
Current flows eastward at about 45°N and splits into the Alaska 
Current and the California Current as it approaches the eastern 
coastal boundary of the Pacific Ocean. This divergence of the 
current varies seasonally and has been demonstrated by satelite- 
tracked drogue buoys. In the winter, more of the buoys’ turn 
southward, and in the summer more of the buoys turn northward 
(Picard and Emery, 1982). 

The AXBT observations from the northwest portion of the 
project area would be expected to show characteristics of the 
Subarctic Water Mass. Those observations from the southeast 
portion should be influenced by the colder offshore waters of the 
California Current. Both the North Pacific Current and the 
California Current have surface speeds of less than one knot. As 
a result, the water in the project area is subject to constant 
climactic conditions and has time to adjust to seasonal variations. 
During the late-fall time period of this experiment, the ocean 
experiences a net loss of heat. This heat loss drives the winter 
convective mixing process and is the major influence on the area's 
vertical thermal structure. 

The physical structure of the ocean depends on both tempera- 
ture and salinity. The principal salinity feature in the north 
Pacific is a permanent halocline found at about 100 m depth. 
Above this level, the water is nearly isohaline and is of relatively 
low salinity. The low salinity water of this region is caused by an 
excess of precipitation over evaporation. Between 100 and 200 m 


depth, the salinity rises about one ppt. 
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The winter cooling season runs from about mid-September to 
about mid-April. Surface cooling and strong wind mixing cause the 
MLD to deepen. If the winter cooling is severe, the residue of 
the decaying seasonal thermocline can descend into the  halocline 
early in the season. The water in the mixed layer above continues 
to cool, and the result is a temperature inversion in the upper 
portion of the halocline. This structure is hydrostatically stable 
because of the relative domination of salinity over temperature in 
controlling the pycnocline (Tully and Giovando, 1963). If the 
winter cooling is not severe, the residue of the mixed layer will 
not reach the halocline during the cooling season. In this case, 
the halocline would remain intact and there would be no temperature 


inversion at depth. 


DEPTH (m) 





STATION NUMBER 
Figure 2.1 Sound veloc y contours aced a 1 m/s showing 
shallow sound channel (after Bem ap et sy A98). 
The existence of a temperature inversion at depth will result in 
the formation of an acoustic shallow sound channel (SSC). The 


sound channel intensity or magnitude is measured in units of sound 


velocity and is the difference between the sound velocity at the 
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axis and that occurring at the lower boundary (Dunlap, 1982). A 
sample set of sound velocity contours is plotted using EOTS data 
for the ASTREX region (Fig. 2.1). The existence of an SSC can 
be seen as the transect extends seaward, from station 8 to 10 at 


about 110 m depth. 


To 


III. MODEL THEORY 


A. OCEAN THERMODYNAMIC MODEL 


The ocean thermodynamic model used in this coupled system is 
the Oceanic Boundary Layer Model (OBL) which was developed by 
Garwood (1977). It is a one dimensional, second order turbulence 
closure, vertically integrated (bulk) model of the ocean surface 
turbulent boundary layer. The model employs the Navier-Stokes 
equation of motion without the geostrophic component, the continuity 
equation for an incompressible fluid, the heat equation from the 
first law of thermodynamics, an analytical equation of state, and a 
two-component vertically integrated turbulent kinetic energy (TKE) 
budget. 

The rate of deepening or shallowing of the mixed layer is 
determined by the dynamics of the entrainment process. The THE 
in the mixed layer erodes the stable water mass found below the 
mixed layer. ‘The entrainment hypothesis is based on the TKE 
budget. Closure for the system of equations is achieved by using 
the bulk buoyancy and momentum equations along with the mean 
turbulent field modeling of the vertically integrated equations for 
the individual TKE components (Rosner, 1985). 

The OBL model differs from earlier models in the following 
ways. First, the amount of wind generated TKE to be used in 
mixing is a function of the ratio of the mixed layer depth (MLD) 
to the Obukhov mixing length. Second,  viscous dissipation is 
dependent on a local Rossby number and separate vertical and hori- 
zontal equations for TKE are used (Garwood, 1977). 

The energy for vertical mixing is provided by both buoyancy 
flux and shear production. However, buoyancy flux is a more 
efficient energy source because of its direct contribution to the 


vertical component of TKE. The buoyancy equation is derived from 


e 


a linearized equation of state along with the conservation of heat 


and salinity equations: 


p = po{1i - a(8 - 80o) - B(S - So) (3.1) 


Buoyancy is defined by: 


b = g(po - p)/po (3.2) 


where: 
9 = temperature 
= salinity 
= gravity 
= density 


= thermal expansion coefficient 


VW R > Q V 


= density coefficient for salinity 

In most open ocean regions in low or mid-latitudes, salinity 
does not have much effect, and temperature is the dominant factor 
contributing to density structure (Miller, 1976). However, by using 
buoyancy instead of temperature, the model is able to allow for 
processes such as precipitation and evaporation which may signifi- 
cantly affect surface buoyancy flux. 


The OBL uses several profiles for initialization: 


1. mixed layer temperature; 
2. mixed layer salinity; and 
3. wind-driven horizontal current. 


The lack of initial salinity and wind-driven current profiles is not a 
significant problem because the model will evolve its own profiles 
within the first diurnal cycle. These profiles are not available 
when the input is from AXBTs. However, the model results are 
not sensitive to initial salinity and current profiles (Davidson and 
Garwood, 1984). 

The time step used in OBL is one hour and there are two 


ways of inputing hourly information into the model. The first 
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method, and the one used in this research, is to supply the heat 
fluxes directly from FNOC. The second method is to calculate 


these heat fluxes each hour using the following boundary condi- 


tions: 

1. net upward turbulent heat flux at the sea surface; 

zi incident solar radiation; 

3.  absorption factor for short wave radiation in the top meter 
of the sea; 

4. wind speed and wind direction; 

di cloud cover; 

6. sea surface temperature (SST); 

R dry bulb air temperature; 

8. dew-point; and 

Ji precipitation. 


& 
Several other physical and model constants are input including: 


1. extinction coefficient; 
2» Coriolis parameter; 
>. critical Richardson number for stability adjustment below the 


mixed layer; 
expansion coefficient for temperature; and 
density coefficient for salt. 
The upward heat flux, Qu, can be calculated by summing the 
turbulent flux of latent heat, Qe, the turbulent flux of sensible 


heat, Qh, and the net back radiation, Qb: 

Qu = Qe + Qh + Qb (359 
The Qe term is calculated using: 

Qe = Cd(.98Es - Ea)U'**”? (3.4) 
where: 

Cd = coefficient of drag 


Es = air vapor pressure of marine air 


Ea = air vapor pressure based on dew point temperature 
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U'!%) = wind speed at 10 meters 


The Qh term is calculated using: 


@h =sed(Ts = Ta)U‘**) (SEO) 
where: 
Ta = air temperature in degrees Kelvin 


Ts = sea surface temperature in degrees Kelvin 
The Qb term is estimated using an empirical equation (Husby and 
Seckel, 1978): 


Ob < 1.14 X 107? (Ts)" (.39 - .5Ea!4)(1 - .6C?) (3.6) 
where: 
C = cloud cover 


The net radiation is not just the simple difference between the 
upward and downward heat fluxes because the shortwave radiation 
is absorbed throughout the mixed” layer (Swaykos, 1989). 
Approximately 50% of the energy is absorbed in the top meter and 
the rest is absorbed below that. Very little energy ‘penetrates 
below the mixed layer because the underlying stable thermocline has 
very low thermal conductivity. The percentage absorbed in the top 
meter (Rf) varies from region to region and is most dependent on 
the amount of suspended particulate matter. Knowledge of this 
absorption factor allows calculation of the net heat flux, Qn, at the 


surface: 
Qn = Qe + (RF)Qs - Qs (321) 


The surface fluxes of momentum and buoyancy can be calcu- 
lated when Qn is known. Mixed layer temperature, salinity, 
velocity and buoyancy fluxes can be computed using four equations 
and four unknowns. These fluxes determine mixed layer shallowing 


or deepening by entrainment. 
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The boundary layer will always deepen if there is a positive 
(upward) surface buoyancy flux. The entrainment heat flux is 
determined from the entrainment model (Garwood, 1977), and it 
determines the rate of mixed layer deepening. This deepening rate 
is used to compute a new temperature profile each time step. The 
boundary layer will shallow if the surface buoyancy flux is negative 
and there is not enough wind mixing and resultant vertical TKE to 
transport heat down to the bottom of the earlier-established mixed 
layer. A new mixed layer will be created at the depth where the 
vertical flux approaches zero. When there is no entrainment and 
the boundary layer shallows, there will be a unique solution for the 


depth of the mixed layer which allows for conservation of heat. 


B. OCEAN ACOUSTIC MODEL 


The ocean acoustic model used in the second part of this 
coupled system is the RAYMODE passive propagation loss model 
developed by G. A. Leibiger at Naval Underwater Systems Center, 
New London Laboratory (Leibiger, 1971). This model is contained in 
the RAYMODE Passive Propagation Loss Computer Program -and the 
version used in this research was written in 1982. The term 
RAYMODE represents a technique of combining aspects of ray theory 
and normal mode theory to cälculate acoustic pressure. Normal 
mode theory allows calculations with more accuracy and detail than 
does ray theory and avoids the ambiguities when rays cross and 
create false caustics. The use of ray theory simplifies some of the 
normal mode calculations because it -uses the usual geometrical inter- 
pretations of ray paths. 

The RAYMODE program uses the sound speed profile (SSP) to 
partition the corresponding wavenumber domain into regions which 
relate to the different propagation paths. The acoustic pressure 
field is calculated by a numerical integration of the normal modes 
for each path. By summing over all of the different propagation 


paths, the total field is obtained. The components of this pressure 
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field are used to form coherent and random phase summations which 
give the total pressure field at each desired range. The difference 
between coherent and random propagation losses is that the coherent 
losses result from the contribution of each ray path assuming the 
actual phase between the rays, while the random losses are summed 
assuming a random phase distribution (Medeiros, 1982). 

RAYMODE has several advantages that make it a good program 
for the at-sea acoustic forecaster who requires sonar performance 
predictions (Yarger, 1982): 

1. This version can be run on small machines such as the HP 

9845 or Tektronix 4051 desktop computers; 
2. It considers and sums together ducted, convergence zone, 
bottom bounce, and other propagation paths; 
It handles a large range of sonar frequencies; 
It accommodates any type of sound velocity profile; and 
It produces valid results for deep or shallow water. 
There are some drawbacks for this 1982 version of RAYMODE 
(T.B. Gabrielson, personal communication, 1985): 

I It does not allow for change of any environmental parame- 
ters along the propagation track even if these changes are 
known; 

2: The leakage from the surface duct is ignored for frequen- 


cies less than 100 Hz; 


3. The integral approximations can be poor for frequencies less 
than 100 Hz; 

4. The smoothing in the incoherent mode is excessive; and 

o. Too few bottom bounce rays are used. 


It is possible that the problem of not being able to input changing 
environmental parameters could be corrected if a range-dependent 
program is developed in the future. 

The original approach (Fisher, 1981) to combining thermody- 
namic and acoustic models used the Fast Asymptotic Coherent 
Transmission Loss model (FACT) for acoustic prediction. Fisher 


stated that the FACT model had several known problems: 
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1. Low frequency problems caused by large scale cancellations 
of caustics due to the inability of ray theory to deal with 
the resulting intensities; 

2. Half channel cases historically presented a problem due to 
surface interference. This problem has been addressed in 
recent updates to the FACT model; however, the the 
existing routine available at the Naval Postgraduate School 
is FACT-9H, and it has not been modified for half channel 
cases; and 

3. The propagation in the surface duct is virtually unaffected 
by wind speed and wave height. For some low frequen- 
cies, unrealistic duct thicknesses are required. 

A comparison between RAYMODE-X, FACT-9H, and two ray 
tracing models was carried out (Hall and Holt, 1981). The basic 
physics behind -the models was discussed and all models were run 
using identical inputs. At low frequencies (below 1 kHz), 
RAYMODE encountered the fewest problems with a surface duct. 
This is because it employs a combination of both normal mode and 
ray tracing inside the duct. The other three ray tracing models 
are forced to perform special calculations for duct propagation. 
The FACT model characterizes gross features of duct propagation 
using the Clay-Tatro equations. These equations use conservation 
of energy along with range-related leakage and scattering losses to 
calculate intensity when both the source and receiver are located in 
the surface duct. For cases when only the source or the receiver 
is in the duct, the intensity is reduced by 10 dB. The use of 
these different equations causes inconsistent results in propagation 
loss for velocity profiles that contain a surface duct. 

The RAYMODE model was chosen for this research because the 
combination of normal mode and ray tracing should, in theory, work 
better for low frequencies and surface duct cases. The OBL 
predicts temperature profiles for the upper layer of the ocean, and 
it has very little effect in the short term on temperature below 150 


m. In coupling the OBL to an acoustic model, the most important 
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results will be those obtained for shallow receivers and  shallow 


sources where the variation of MLD has the greatest effect. 
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IV. MODEL RESULTS 


A. THERMODYNAMIC MODEL RESULTS 
1. General 


The thermodynamic model (OBL) was initialized using the 
AXBT data acquired 15 November 1980 during the  STREX field 
phase. The stations used were 2, 3, 4, 5, 6, 7, and 8 (see Fig. 
1.1) because they were located on the center track, and that track 
was the outbound leg on all P-3 flights. No data were available 
for station 1 as the 15 November AXBT did not deploy properly. 
The OBL was run for 20 days and the model output was compared 
with final AXBT data acquired 5 December 1980 at these same 
stations. 

The model used atmospheric forcing which had been 
obtained from FNOC and placed in mass storage files. A correction 
factor was applied to the supplied heat fluxes to account for Known _ 
bias in the forcing- (Gallacher, 1979). The program interpolated to 
the closest two degrees of latitude and five degrees of longitude to 
obtain wind speed, solar radiation, and heat flux. An example of 
this atmospheric forcing which was interpolated for station 2 is 
plotted for the 20 days (Fig. 4.1). The model used a time step 
of one hour and integrated for 20 days. Each station will be 
discussed independently and then the stations will be considered 


together for a two dimensional analysis. 
2. Station 2 


The OBL was initialized at station 2 with AXBT data that 
indicated a MLD of 32 m and a surface temperature of 13. 8 CH CEE 
4.2): After integrating for 20 days, it predicted that due to 
increased TKE, the MLD would deepen to 58 m, and the surface 


temperature would cool 2.6°C to 11.2°C, The final AXBT RL TE 
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Figure 4.1 Interpolated atmospheric forcing for Station 2. 


days later showed a MLD of 07 m. However, the surface tempera- 
buče cooled only 0.79€ to 13.1°C. Of the seven stations studied, 
this was by far the smallest amount of surface cooling and was the 
only station to show a local temperature maximum at 125 m. A 


local temperature maximum at this depth could be caused by 
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Figure 4.2 Comparison between AXBT and OBL at Station 2. 


horizontal advection or by a deep mixed layer earlier in the year. 
Since there was no evidence of this feature in the AXBT taken 20 
days before, the results suggest that this station was either 


affected by horizontal advection or influenced by a mesoscale feature 


which had been missed (navigational error) by the previous AXBT 
observation. The one dimensional OBL initialized with AXBT data is 


not able to resolve either of these situations. 


INITIAL AXBT FINAL AXBI 


6 7 8 9 10 LI 12 13 14 15 16 6 7 8 9 10 11 12 13 14 15 16 
0 0 


6 7 8 9 10 11 12 13 14 15 18 6 7 8 9 10 11 12 13 14 15 18 
TEMPERATURE (*C) TEMPERATURE (°C) 


INITIAL MODEL FINAL MODEL 


6 7 8 9 10 11 12 13 14 15 18 6 7 8 9 10 11 12 13 14 15 16 
0 


6 7 8 9 10 11 12 13 14 15 18 6 7 8 9 10 11 12 13 14 15 16 
TEMPERATURE (°C) TEMPERATURE (°C) 





Figure 4.3 Comparison between AXBT and OBL at Station 3. 
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3: Station 3 


The OBL at station 3 was initialized using AXBT data 
which had a MLD of 66 m and a surface temperature of 15.9°C 
(Fig. 4.3). The model was integrated for 20 days and predicted a 
MLD of 77 m with a surface temperature of 14.7°C. The final 
AXBT 20 days later showed a surface temperature of 14.4°C and a 
MLD that only deepened very slightly to 68 m. At this station, it 
is possible that the model calculated more TKE than actually 
existed. The shape of the temperature profile below the layer 
shows a very good correlation between the model output and the 
final AXBT information, with the exception that the model output 
tends to be smoother. This smoothness is to be expected as the 
model employs a weak below-layer diffusion and does not account 
for the small-scale horizontal currents and internal wave motion that 


influence the actual temperature profile. 
4. Station 4 


The OBL at station 4 was initialized with AXBT data that 
had a MLD of 52 m and a surface temperature of 14.9 CA 
4.4). The model was integrated for 20 days and predicted a MLD 
of 67 m with a surface temperature of 13.5°C. The final AXBT 
data from 20 days later, showed a surface temperature of 13.8°C 
down to the MLD of 72 m. The profile shapes of the model output 
and the final AXBT information look very similar with the major 
difference again being that the model profile is much smoother. An 
interesting feature of the. initial AXBT conditions was that the 
cooling of the top part of the mixed layer appeared to have already 
started. The cooler water at the surface at this station is 
unstable because the salinity is well mixed in the mixed layer. 
This initial trace clearly shows the cooling source of the TKE for 


mixed layer deepening. 
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Figure 4.4 Comparison between AXBT and OBL at Station 4. 


5: Station 5 


The OBL at station 5 was initialized with AXBT data that 
showed a surface temperature of 14.4°C down to a MLD of 42 m 


(Fig. 4.5). The model was integrated for 20 days and predicted 
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Figure 4.5 Comparison between AXBT and OBL at Station 5. 


a surface temperature of 12.5°C down 


to a MLD of 64 m., The 
final AXBT taken 20 days later 


showed a surface temperature of 
13.0°C down to a MEDE 968 vin. This AXBT taken 5 December 


was never properly digitized and was not part of the STREX data. 
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This trace was recovered from the original aircraft files ang read 
using a metric AXBT grid overlay. The profile could only be 
resolved to about 0.2°C and to about 3 m in depth. Even with 
this resolution, the AXBT data confirmed that the model deepened 
to approximately the correct depth, but showed that the model 
surface temperature was too cold by 0.5°C. This could be due to 
inaccuracies in the heat fluxes provided by FNOC or the flux field 
not being ,resolved closely enough for the particular station. The 
difference could also be due to surface currents which would result 


in horizontal advection and a relative warming of the mixed layer. 
6. Station 6 


The OBL at station 6 was initialized with AXBT data that 
showed a surface temperature of 14.3°C down to a MLD of 63 m 
(Fig. 4.6). The model was integrated for 20 days and predicted 
a surface temperature of 12.7°C down to a MLD of 78 m. The 
final AXBT taken 20 days later showed a surface temperature of 
12.6°C down to a MLD of 77 m. The OBL predictions made at this 
station were very accurate. Almost every point on the profile 
down to 200 m was verified within 01°C: Perhaps this station was 
not affected by any significant mesoscale activity during the 20 


days. 
7. Station 7 


The OBL at station 7 was initialized with AXBT data that 
showed a surface temperature of 13.4°C down to a MLD of 49 m 
(Fig. 4.7). The model was integrated for 20 days and predicted 
a surface temperature of 11.3°C down to a MLD of 72 m. The 
final AXBT taken 20 days later showed a surface temperature of 
11.8%C down to a MLD of 72 m. The OBL again predicted the new 
MLD within a meter but showed 0.5°C too much surface cooling. 
The final AXBT profile below the layer looked very similar to the 
OBL profile with the exception of a small region between 150 and 
170 m. The AXBT profile showed a very slight positive gradient 
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Figure 4.6 Comparison between AXBT and OBL at Station 6. 


at 200 m depth that did not exist 20 days before and could not be 
developed by the model. This must be due to horizontal advection, 
internal vertical motion of the water column, or movement of a 


sub-surface mesoscale feature with respect to the station position. 
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Figure 4.7 Comparison between AXBT and OBL at Station 7. 


Another possibility is that the AXBT observations 20 days apart 


were not in the same position. 
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Figure 4.8 Comparison between AXBT and OBL at Station 8. 


B. Station 8 


The OBL at station 8 was initialized with AXBT data that 
showed a surface temperature of 13.3°C down to a MLD of 57 m 


(Fig. 4.8) The model was integrated for 20 days and predicted 
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a surface temperature of 11.2%C down to a MLD of 83 m. The 
final AXBT taken 20 days later showed a surface temperature of 
11.2°C down to a MLD of 77 m. The OBL and AXBT profiles 
appeared very similar at this station with the only exception being 
the smoothness of the OBL profile. The predicted MLD was 6 m 
too deep, but the temperatures in the upper 200 m were almost all 
verified within 0.2°C. Both profiles had isothermal regions between 
150 and 180 m, and that feature was also apparent in the initial 
AXBT. 


9. Two Dimensional Analysis 


The main reason that this ASTREX AXBT information was 
selected for analysis was that the one dimensional model could be 
run at several equally spaced points along a track. This enabled 
an easy transition from one dimension to two dimensions. The 
initial AXBT information was plotted from station 2 to station 8 and 
contoured for temperature down to 300 m (Fig. 4.9). These 
contours provided a view of the synoptic features which were not 
readily apparent from the individual AXBT traces. Each station is 
separated by about 50 km and the track heads northwestward from 
the Californian coastline. 

le oo ima ostal (etui renuts are apparent as the surface 
waters just off the coast are initially cool and become warmer 
between stations 2 and 3. The MLD is deepest (67 m) at station 3 
and every isotherm in the transect is at its deepest point here. 

After the OBL has integrated for 20 days (Fig. 4.10), the 
predicted temperature contours are plotted for station 2 to station 
8. The MLD has deepened at all stations, and the temperature in 
the layer has cooled an average of 1.9°C. The surface cooling was 
the greatest where the mixed layer was the shallowest, because the 
colder water below was entrained into the mixed layer as the TKE 
increased due to surface cooling and wind stress. 

The final AXBT contours are plotted for the same stations 


(Fig. 4.11). The MLD deepened as predicted and is well defined 
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Figure 4.9 Initial AXBT contours x temperature versus depth 
15 November 1980. 


in this transect. The surface layer cooled an average of 1.4°C, or 


0.5°C less than wás predicted This difference in surface 


13.0 
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Figure 4.10 Predicted contours oj temperature versus depth 
5 December 1980. 


temperature was seen in the analysis at the individual stations and 


may have heen as a result of the heat fluxes supplied by FNOC. 
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Figure 4.11 Final AXBT contours of temperature versus depth 
December 1980. 


Although these fluxes included a correction factor (Gallacher, 1979), 


it.is possible that a more accurate flux field could be calculated for 


this portion of the Pacific. 
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The warmest mixed layer water is found at station 3. 
The warm synoptic feature predicted for this station down to 250 m 
was verified, but only down to 150 m. A temperature inversion 
can be seen at station 4, and the final AXBT data contain other 
inversions but they are not large and are not resolved by the 
contours of Fig. 4.11. The 10°C isotherm in both transects is 
found at the bottom of the steep temperature gradient found imme- 
diately below the mixed layer. It appears in almost the same place 
in both the predicted and observed  transects. Both sets of 
contours show a convergence of isotherms just below the MLD as 
the transect extends seaward in the northwestward direction. 
Overall, these transects demonstrate that the OBL, if initialized 
properly and forced by the atmosphere,  posseses a significant ` 
predictive capability. 

& 

B. ACOUSTIC MODEL RESULTS 


e General 


The RAYMODE model was run three times at each of the 
seven stations. The first run was made using the initial AXBT 
information taken 15 November 1980. The second run was made 
using the thermal structure calculated for 5 December 1980 after the 
model was integrated for 20 days. A final run was made using the 
AXBT information available for these stations on 5 December. No 
changes were made to any of the RAYMODE parameters or inputs 
except for the temperature profile. The RAYMODE program uses 
salinity and depth to calculate a velocity for each point that a 
temperature is input. In all three runs, a standard historical 
salinity profile was used at all stations and was not changed for 
the 20 days. Neglect of salinity variations has a negligible effect 
on sound speed calculations in this region of the Pacific Ocean. 
The major effect of salinity variations is that changes in the 


salinity profile permit temperature inversions. 
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The user has the choice whether to use coherent or 
random propagation loss information. The coherent output is theo- 


retical and shows large fluctuations every tenth of a kilometer 


which is not realistic for the operational user. In this research, 
the random output is used for both the MDR and the CZR calcula- 
tions. The random output has been  smoothed for the user. 


However, the amount of smoothing carried out is high and this is a 
recognized problem with the RAYMODE program. 

The acoustic results were based on one frequency (300 
Hz), one wind speed (5 kts), one receiver depth (20 m), and one 
source depth (60 m). It is recognized that all of these variables 
could be changed and studied with each of the temperature inputs, 
but a sensitivity analysis of the RAYMODE model is beyond the 
scope of this research. The shallow source and receiver were 
selected because the OBL models the upper 200 m of the ocean and 
finds the most variability in the upper 80 m. The 20 m receiver 
depth was selected because it was very close to the standard 
shallow sonobuoy depth of 60 feet. The 60 m source depth was 
selected because it is of interest to model a scenario where the 
submarine is acoustically hiding just below the mixed layer. The 
frequency of 300 Hz was chosen because it was high enough to 
avoid the low frequency problems experienced by the RAYMODE 
model and low enough to allow good passive acoustic information to 
be obtained. The arbitrary wind speed of 5 knots was selected 
because that low wind speed causes very little interference in the 
surface duct and at 5 knots, the RAYMODE Program will not 
experience any problems as the wind speed approaches Zero. 

Two measures of acoustic performance were used in this 
evaluation of acoustic output. The first was median detection range 
(MDR). This was calculated from the propagation loss curve at 
figures-of-merit (FOM) OL (oR 80, and 85 decibels (dB). 
Operationally, a user in the fleet enters a propagation loss plot 
with an FOM and extracts the MDR from the intersection of the 


curve and the FOM. Three values of FOM were used to lend 
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TABLE I 
RAYMODE Output for MDR and CZR in Kilometers 


Station 2 
FOM (dB) INITIAL MODEL AXBT 
MDR CZR MDR CZR MDR CZR 
75 2.0 - 205 - 28 - 
80 3.3 - 3.4 32953 3.8 35.8 
85 5.1 35.7 AT 30.8 54 34.6 
Station 3 
FOM (dB) INITIAL MODEL AXBT 
MDR CUR MDR CZR MDR OZR 
75 25 - 2.4 - 2.6 - 
80 oe - 33 38.6 Sed Set 
85 45 39.1 4.3 Jia 4.6 36.2 
Station 4 
FOM (dB) INITIAL MODEL AXBT 
MDR CZR MDR SR MDR SR 
75 2.9 - 2.6 - Zak - 
80 > 37.9 3.5 35.6 3.6 36.1 
85 4.5 36.9 4.9 34.4 Sal 2I 72 
Station 5 
FOM (dB) INITIAL : MODEL AXBT 
MDR CZR MDR OZR MDR GER 
19 3 = DES = 2.6 - 
80 pr - 3.4 37.9 3.5 36.1 
85 4.3 21.3 4.8 33.6 4.8 34.3 
Station 6 
FOM (dB) INITIAL MODEL AXBT 
MDR GIR MDR CZW MDR CZR 
75 25 - DINO, - 2.4 - 
80 3.4 37.9 320 35.4 3.3 35.6 
85 4.6 36.1 4 2 3016 4.3 Sout 
Station 7 
FOM (dB) INITIAL MODEL AXBT 
MDR CZR MDR GER. MDR CAR 
75 AO > 2.4 = De - 
80 Sal = 22 Soc Siren, zana 
85 a 36.0 4.4 3159 4.4 SA 
Station 8 
FOM (dB) INITIAL MODEL AXBT 
MDR CZR MDR G7 i. MDR GZR 
15 205 = 2.4 = 24 = 
80 Dd 35 SS 34.3 Se oa. 2 
85 4.6 Do 2 0: 4.4 32.6 4.3 33.0 
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insight to trend analysis. The MDR is an important acoustic model 
output because it is used by ASW fleet units such as patrol 
aircraft, helicopters, and towed array ships for sonobuoy spacing 
and for tactical planning. 

The second measure of acoustic effectiveness used was 
convergence zone range (CZR). A fleet user obtains CZR from the 
same propagation loss plot and continues further along the FOM line: 
until there is an intersection again with the propagation loss curve. 
The CZR is also an important acoustic model output because knowl- 
edge of the expected CZR allows for faster resolution between 
direct path and  convergence zone paths from passive acoustic 
receivers. The lack of, or change in expected CZR can have 
important tactical implications on the deployment of sonobuoys and 
ASW aircraft. The bottom depth along the track ranged from 3000 
to 4000 m as the stations extended seaward with the deep sound 
channel between 250 and 300 m. There was more than  sufficient 
depth excess at all stations to insure CZ propagation. There is 
not a convergence zone ‚in all cases, and in this particular anal- 
ysis, the 75 dB FOM did not give a CZR at any of the stations. 
The MDR and CZR outputs have been tabulated (Table I) for the 
RAYMODE (runs. | 


da Median Detection Range 


The MDR is plotted at all stations for the three different 
FOM cases studied (Fig. 4.12). As would be expected, the MDR 
increases in both range and variability with increasing FOM. The 
amount of MDR variability between the three runs at each station is 
low, but it is clear from the plot that the model predicts the MDR 
well. 

At an FOM of 75 dB, the effect of a shallow MLD and a 
steep gradient below the layer is clearly seen at station 2. After 
the MLD has deepened, the MDR is more consistent with the results 
at the other stations. With an FOM of 80 dB, there is more varia- 


tion over the 20 day. time period, and again the model predicts 
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STATION NUMBER 





Figure 4.12 RAYMODE MDR for FOM of 75, 80, and 85 dB.. 


MDR very well. The case of an FOM of 85 dB shows the most 
variability over the 20 days, and the model predicts the MDR well 
except at station 2. This may be due to the fact that the OBL 
was off by nearly two degrees in the layer temperature and missed 


the below layer temperature maximum. 
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STATION NUMBER 





Figure 4.13 RAYMODE CZR for FOM of 80 and 85 dB.. 


s Convergence Zone Range 


The CZR is plotted at all .the stations for the two 
different FOM cases studied. There was no convergence zone (CZ) 
below 77 dB at any of the stations in any run using the random 
propagation loss output. There is clearly much more variability in 
CZR than in MDR. This variability is to be expected because small 
changes in surface temperature are amplified by the deep conver- 
gence zone rays. 

At an FOM of 80 dB, the initial conditions. showed A 
only at stations Io ande: The model output predicted a CZ at 
all the stations as did the computations using the AXBT data. The 
variability associated with the model run at station 2 and the warm 


synoptic feature at station 3 is evident. At an FOM of Bozen 
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this variability at station 2 is also evident as is a clearly defined 
maximum CZR at station 3. The range where the CZ is first 
detected is well predicted by the model and, at all stations, is less 
than the range given by the initial conditions. This is due to 
surface cooling where both the source and receiver are located, and 


the fact that the deep ocean characteristics remain unchanged. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 


The coupled mixed layer-acoustic model was subject to analysis 
at a line of stations in the northeast Pacific Ocean. Since no 
real-time acoustic data were available for verification, the OBL ther- 
modynamic model output was compared to final AXBT data after both 
cases were run through the RAYMODE acoustic model. Thermal and 
acoustic analysis of the data led to several important conclusions: 

1. The OBL model, after integrating for 20 days, accurately 

predicts the magnitude of the deepening of the mixed layer; 

Die The coupled model accurately predicts trends in MDR and 
CZR by determining MLD and surface temperature using the 
physics of entrainment; 

3, The variability of MDR during the 20 days was very low 
and demonstrates how slowly the thermal structure changes 
during the late fall period and how important it is to 
initialize acoustic models with actual data from the exact 
location being studied; and 

4. The OBL, and hence the coupled system, produces better 
results when the station at which it is run is not influ- 
enced by near-shore currents or by horizontal advection of 


any kind. 


B.  RECOMMENDATIONS 


This coupled model could be used operationally in the fleet and 
would be an excellent aid to the on-board tactical commander. The 
situation that presently exists in the fleet calls for the user to 
rely on historical BT profiles generated for arbitrary grid squares 
when a current profile is not available. Even if a fairly recent 


profile is available, historical data bases currently in use at FNOC 
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such as the Expanded Ocean Thermal Structure (EOTS) model 
discard actual BT information after an arbitrary time. Many situ- 
ations arise where a fairly recent profile may be obtained, but no 
method is available to advantageously apply that information because 
surface conditions such as temperature and MLD may have changed. 
The following recommendations are made to enable the coupled model 
to be used operationally: 

1. The OBL and RAYMODE programs could be combined so that 
both models can be run using the output of one as the 
input for the other; 

2, This coupled program should be adapted for use on the 
small computers presently employed for acoustic prediction at 
sea; 

3. A method should be found to allow the fleet user easy 
access to the atmospheric forcing files‘ available from ENOC. 
and 

4. Further acoustic analysis of the coupled model output should 
be carried out. The OBL has been throughly tested and 
is fairly reliable. However, its interaction with an acoustic 
model needs to be studied for different acoustic and ther- 
modynamic scenarios. 

This coupled model has much potential for the tactical 
commander in the fleet. It offers a predictive capability that does 
not exist now on shore or at sea. Adapted for use at sea, this 
system would produce reliable short term acoustic forecasts in a 
stand-alone mode. The task force then would be able to make use 
of AXBT information for a specific area even if that information was 


several days or weeks old. 
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